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Effects of Resonance in Quasiperiodic Oscillators of Neutron Star 

Binaries 

Lev Titarchuk 1,2 
ABSTRACT 

Using a large quantity of Rossi X-ray Timing Explorer data presented in the 
literature I offer a detailed investigation into the accuracy of quasiperiodic oscil- 
lations (QPO) frequency determination. The QPO phenomenon seen in X-ray 
binaries is possibly a result of the resonance of the intrinsic (eigen) oscillations 
and harmonic driving forces of the system. I show that the resonances, in the 
presence of the damping of oscillations, occur at the frequencies which are sys- 
tematically and randomly shifted with respect to the eigenfrequencies of the 
system. The shift value strongly depends on the damping rate which is measured 
by the halfwidth of the QPO feature. Taking into account this effect I analyze 
the QPO data for four Z-sources: Sco X-l, GX 340+0, GX 5-1, GX 17+2 and 
two atoll sources: 4U 1728-34, 4U 0614+09. The transition layer model (TLM) 
predicts the existence of the invariant quantity: 5, an inclination angle of the 
magnetospheric axis with respect to the normal to the disk. I calculate 5 and 
the error bars of 5 using the resonance shift and I find that the inferred 5— values 
are consistent with constants for these four Z-sources, where horizontal branch 
oscillation and kilohertz frequencies have been detected and correctly identified. 
It is shown that the inferred 5 are in the range between 5.5 and 6.5 degrees. I 
conclude that the TLM seems to be compatible with data. 

Subject headings: Accretion, accretion disks — stars: individual (Sco X-l, GX 
340+0, GX 5-1, GX 17+2, 4U 1728-34, 4U 0614+09) — stars: neutron 
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1. Introduction 

Kilohertz quasi-periodic oscillations (QPOs) have been discovered by the Rossi X-ray 
Timing Explorer (RXTE) in a number of low mass X-ray binaries (Strohmayer et al. 1996, 
van der Klis et al. 1996). The presence of two observed peaks with frequencies v\ and 
z/ 2 in the upper part of the power spectrum became a natural starting point in modeling 
the phenomena. Attempts have been made to relate v\ and z/ 2 and the peak separation 
Az/ = z/ 2 — v\ with the neutron star (NS) spin. In the sonic point beat frequency (SPBF) 
model by Miller, Lamb & Psaltis (1998) the kHz peak separation Az/ is considered to be close 
to the NS spin frequency and thus Az/ is predicted to be constant [see also the updated version 
of the SPBF (Lamb & Miller 2001) in which Az/ is allowed to be variable]. Observations of 
kHz QPOs in a number of binaries (Sco X-l, 4U 1728-34, 4U 1608-52, 4U 1702-429 and etc) 
show that the peak separation decreases systematically when kHz frequencies increases (see 
a review by van der Klis 2000). 

There are two other models in the literature which infer the relations between v 1 , z/ 2 
and vhbo- The relativistic precession (RP) model involves high speed particle motion in 
strong gravitational fields, leading to oscillations of the particle orbits. Stella et al. (1999) 
studied precession of the particle orbit under influence of a strong gravity due to General 
Relativity (GR) effects. The correlation between high frequency (lower kHz frequency) and 
low frequency (broad noise component) QPOs previously found by Psaltis, Belloni & van der 
Klis (1999) for black hole (BH) and neutron star (NS) systems has been recently extended 
over two orders of magnitude by Mauche (2002) to white dwarf (WD) binaries. With the 
assumption that the same mechanism produces the QPO in WD, NS and BH binaries, 
Mauche argues that the data exclude RP, SPBF models (as well as any model requiring 
either the presence or absence of a stellar surface or a strong magnetic field). 

The transition layer model (TLM) was introduced by Titarchuk, Lapidus & Muslimov 
(1998), hereafter TLM98, to explain the dynamical adjustment of a Keplerian disk to the 
innermost sub-Keplerian boundary conditions (e.g. at the NS surface). The low branch 
frequency i>l, the Kepler frequency z/ K and the hybrid frequency u h) as they are introduced by 
OT99 are eigenfrequencies of the oscillator. In the TLM framework the observed horizontal 
branch oscillation (HBO) frequencies vhbo, v \ and z/ 2 can be interpreted as the resonance 
frequencies near u Ll z/ K and u h which are broadened as a result of the (radiative) damping 
in the oscillator (see TLM98, Eq. 15). Furthermore, the resonance frequencies vhbo-, v \, v 2 
are shifted with respect to the eigenfrequencies z/£, z/k, Vh- The frequency shift and random 
errors of the eigenfrequencies depend on the damping rate of oscillations A (see details in §3). 
One should keep in mind that the systematic and random errors in the centroid frequency 
determination due to this resonance shift can be a factor of a few larger than the statistical 
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error in the determination of the centroid frequency. It is worth noting that an interpretation 
of the observed QPO phenomena in the framework of any oscillatory (or wave propagation) 
model requires the aforementioned corrections due to resonance in the presence of damping. 
For the observed Q- values (Q = v/2X) of a few, the damped intrinsic harmonic oscillations 
decay very quickly (less than 10~ 2 s) and thus only the forced oscillations can exist for time 
periods of order 10 2 ~ 3 seconds and longer (for details see §3.1). There remains the question 
of what kind of the observational arguments can support the resonance mechanism for QPO 
phenomena. I address to this issue in §4. The goal of this Letter is to demonstrate how the 
resonance shift corrections are important for the accurate interpretation of the QPO data. 
The Letter reports the results of the QPO data interpretation from four Z-sources: Sco X-l, 
GX 340+0, GX 5-1, GX 17+2 and two atoll sources: 4U 1728-34, 4U 0614+09 and two atoll 
sources, collected by RXTE during 6 years of observations. In §2, I present the origins of the 
RXTE QPO data used in the present investigation. In §3 I provide details of the resonance 
effect on the eigenfrequency restoration using the observed QPO frequencies and in §4 I offer 
comparisons of the predictions of the TLM model with the RXTE observations. Summary 
and conclusions are drawn in §5. 



2. Source and Data Selection 

I used the QPO data from Sco X-l [van der Klis et al. (1996), (1997), Bradshaw et al. 
(2002)]; from GX 340+0 (Jonker et al. 2000a); from GX 5-1 (Jonker et al. 2002); from GX 
17+2 (Homan et al. 2002); from 4U 0614+09 (van Straaten et al. 2000, 2002) and from 4U 
1728-34 (van Straaten et al. 2002). 



3. Resonance effects in weakly nonlinear oscillators: relevance to observed 

QPO phenomenon in X-ray binaries 

The observational appearance of asymmetric Lorentzian QPO features, HBO and low 
frequency harmonics (e.g. van der Klis 2000) is the first indication of the classical resonance 
phenomena along with a combination frequency effect well established in weakly nonlinear 
oscillating systems (e.g. Landau & Lifshitz 1965, hereafter LL). Similar effects are described 
by radiation theory in the context of atom excitation (as the classical damped harmonic 
oscillator) which is driven by a harmonic wave (see e.g. Mihalas 1970). 

OT99 formulate the QPO problem in the framework of an oscillator in the rotational 
frame of reference. Taking the magnetospheric rotation with angular velocity fl = 27rz/ mg 
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(not perpendicular to the plane of the equatorial disk) the small amplitude oscillations of the 
fluid element thrown into the magnetosphere are described by a three dimensional oscillator 
(see OT99, Eqs 2-4). When the angle 5 , between f2 and the vector normal to plane of 
radial oscillations is small then three dimensional intrinsic oscillations are split into two in- 
dependent eigenmodes: the radial mode with the hybrid eigenfrequency Vh = (v\ + Av^ ng ) 1 / 2 
and the vertical mode with the low branch eigenfrequency v L = 2b / mg (b / K /b / h ) sin<5. Thus 
one can consider two independent harmonic oscillators under influence of a damping force 
probably the radiation drag force (proportional to the velocity of displacement) as intro- 
duced by TLM98. The restoring force in these oscillators is only in the first approximation 
proportional to the dispacement with the proportionality coefficient equal to the square of 
the eigenfrequency. The interplay between the gravitational, Coriolis forces and pressure 
gradient maintains the oscillations with frequencies close to Vh and vl but not exactly at 
these frequencies. In the next approximation one needs to introduce the nonlinear terms for 
the restoring force (see LL and Titarchuk 2002). Thus the small amplitude oscillations of 
the fluid element in the magnetosphere can be described by two equations 

x + 2\ii + u\x = Aie tu)t + ai^x 2 + 2a^ 2 xz + ai^z 2 , (1) 

z + 2\ 2 z + u\z = A 2 e iuJt + a 2>1 x 2 + 2a 2>2 xz + a 2>3 z 2 , (2) 

where x and z are the radial and vertical components of the displacement vector respectively. 
Here we present equations (2-4) from OT99 (written in Cartesian coordinates) using the 
eigenvector basis (radial and vertical eigenvectors). Ai and A 2 are the damping rates for x 
and z components; u: h = 2nu h , uj l — 2nu L are angular hybrid and low branch frequencies 
respectively (OT99), a; is a frequency of the harmonic driving force and Ai, A 2 are amplitudes 
of x and z components of the driving force, ak,i, (k — 1, 2, I — 1, 2, 3) are coefficients at the 
weakly nonlinear terms of the restoring force. 



3.1. Harmonic oscillations 

The linear terms of the solution of equations (1-2) are 

x = C ltl f-(t,uj h , Ai) + C 1>2 f+(t,u h , Ai) + Bi(u L ) exp (iut), (3) 

z = C 2 ,if-(t,u; L , A 2 ) + C 2t2 f + (t,u L , A 2 ) + B 2 (u h ) exp (iut), (4) 

where f±(t,u*, A) = exp {-[At ± it(u 2 - A 2 ) 1/2 ]}, and B k (u*) = A k /[u 2 - ui 2 + 2i\ k u\. The 
first two terms in Eqs. (3-4) represent the damped harmonic oscillations (it is assumed that 
tf*/Afc > 1). They decay very rapidly within a time interval 1/A& which is of order of < 10~ 2 
s. In fact, a typical measured A^ = \k/2n are always higher than a few hertz (see details 
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below). Thus only the forced oscillations can exist for time periods of order 10 2 3 seconds 
and longer. The absolute values of the forced oscillation amplitudes 



>k 



(2 7 r) 2 [(^-^) 2 + 4A^ 2 ] 1 /2' 



k = 1,2 (5) 



(where B k = b k e ie , = u^, for k — 1, 2 respectively and v = uj/2tt) has a maximum 
at v r k = {y1 k — 2A|) 1//2 . Below we omit the subscript k for simplicity of the presentation 
having in mind the particular case for the eigenfrequency u . Hence the frequency shift of 
the resonance frequency v T with respect to the eigenfrequencies u because of damping is 

r / =K 2 -2A 2 ) 1 / 2 -z/ . (6) 

We find that the maximum of the main resonance amplitude for the linear oscillations is 
not precisely at the eigenfrequency uq, but rather it is shifted to the frequency v r = uq + rj. 
For a linear oscillator v r depends on the damping rate A. From Eq. (6) it follows that 
7] ps —X 2 /u = —X/2Q when Q = u /2\ ^> 1. It is easy to show (see also LL for details) that 
the half- width of the resonance Auhw ~ A for Q ^> 1. Taking a differential from the left 
and right hand sides of the formula for v r one can derive a formula for the random error of 
the eigenfrequency determination Az/ due to the resonance effect: 

Az/ = UrAUr/uQ + 2AAA/z/q- (7) 



3.2. Anharmonic oscillations and p:q resonance in the QPO sources 

The resonance in weakly nonlinear systems occurs at the eigenfrequencies of the system 
uq when the frequency of the driving force is 7 ~ p^o/q and p, q are integers (LL). The 
observed harmonics and subharmonics of HBO frequencies (Jonker et al. 2000a, 2000b 
2002) are probably results of the resonance effect in the weakly nonlinear system. The main 
resonance peak power (for p = 1 and q — 1) is the strongest among all the harmonics because 
the peak power diminishes very quickly with the increase of p and q (LL). The maximum 
of the main resonance amplitude for the nonlinear oscillation is also not precisely at the 
eigenfrequency u . In the nonlinear case the systematic and random shifts are (in addition 
to the damping) affected by the amplitude of the driving force and the relative weights 
(coefficients) a k ,i of the nonlinear terms of equations (1-2), (see LL for details). 

Recently, Abramovicz et al. (2002) drew attention of the community to the possible 
presense of z/i/z/ 2 = 2 : 3 resonance for two kHz frequencies. A similar effect was also 
pointed out earlier by Abramovicz & Kluzniak 2001; Remillard et al. (2002) for BH sources. 
Possibly, it is not by chance these kHz frequencies are detected in the observations because 
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in anharmonic oscillator v\ and v 2 would excite each other (if they are related through 
this 2 : 3 ratio, see above). The oscillations with the frequency v 1 (7) would be a driving 
force for the oscillations with the eigenfrequency z/ 2 (uq) and vice versa. Abramovicz et al. 
(2002) argue that it is hard to escape the conclusion that a resonance is responsible for 
distribution of frequency ratios for Sco X-l and other kHz QPO sources. But a 2:3 ratio 
is a natural consequence of the TL model where z/ 2 ps and v\ ~ uk. In fact, u^/uh = 
[1 + {2v mg /v- K )' 2 \~ 1 / 2 . The magnetospheric rotational frequency v mg should very close to the 
NS spin frequency u NS = 300 — 400 Hz (at least within 15 %). But 2v mg ~ 2u NS are very 
close the observed lower kHz frequency v\ z/ K ) and thus the ratio of u K /u h ps 2 -1 / 2 ps 2 : 3. 

The hybrid frequency 1^, introduced and derived in OT99, is a common feature of 
oscillations in a rotational frame of reference under the influence of a gravitational force. 
When a gas (or fluid) has a nonuniform distribution of mass density, it may exhibit the 
Rayleigh- Taylor instability. Hide (1956) and Chandrasekhar (1961) presented the dispersion 
relation for a fluid rotating with angular velocity £1 (u r = Q./2-k) v^\Y — 4z/ 2 /V^] 1//2 = uo 2 . 
It is easy to show that uq is very close to the Kepler frequency and thus one can find that 
Uh, as a root of the Hide-Chandrasekhar relation, is very close to that obtained by OT99. 
The detailed analysis of the Rayleigh- Taylor instability and oscillations for various density 
distributions, configurations and boundary conditions is given in Titarchuk (2002). 



4. Results and Discussion 

It was shown in the previous section that observed QPO frequencies vhbo, v \ and v 2 
could be treated as the resonance frequencies related to the eigenfrequencies of the linear 
damped oscillators. One can restore a particular eigenfrequency z/ , its error bar Az/ using 
formulas (6)-(7) and values of v r , Au r , a QPO halfwidth, Avhw — A. Thus for each set 
of vrbo ± Auhbo, v\ ± Av>i and v 2 ± Az/ 2 we obtain vl ± Aul, ± Auk and Vh ± Avh 
as restored eigenfrequencies. According to formula (10) in OT99 v mg = {y\ — z/^) 1 / 2 /2 is 
the magnetospheric rotational frequency, which, for highly conductive plasma in the case 
of corotation, should be approximately constant v mg = v mg $. In Figures 1 I present v mg 
versus vk for four Z-sources: Sco X-l, GX 340+0, GX 5-1, GX 17+2. Indeed one may notice 
that v mg = Vmgfl = const is consistent for two Z-sources, Sco X-l and GX 17+2. If the 
magnetosphere corotates with the neutron star (solid-body) rotation, then this procedure 
of v m g— extraction determines the spin rotation of the star. For the other Z— sources, GX 
340+0, GX 5-1 one can find the features of the differential rotation which depend on the 
magnetic field B (see OT99). These features of the differential rotation can also be identified 
for atoll sources 4U 1728-34 and 4U 0614+09. The "real" errors bars Av mg = {v h Av h + 



-7- 



z/kAi/k) I^Vmg of z/ m3 become larger due to these corrections and consequently it is very 
difficult to identify significant deviations of the rotation of the magnetosphere from the 
solid-body rotation (at least for two analyzed Z-sources). The 5— value is calculated using 
formula (1) in TO01, 

5 = arcsin [{v 2 h - 4)~ V W^)]- (8) 
The errors of 5 is evaluated using the differential 

A5 = 5[du L /u L + dv h jv h + du K /u K + (v h dv h + v K dv K ) / Av^]. (9) 

In Figure 2 the dependence of the angle 5 between a normal to the disk and a magneto- 
spheric axis on z/k is shown for the four Z-sources The model invariants should be the same, 
independent of the eigenfrequency branch. Indeed, for each of these Z-sources the 5— values 
show little variation with z/k, Vh, an d vl- within the inferred error bars. The reduced xfed 
for the best-fit 5 is always very close to 1. Jonker et al (2002) presented the analysis of the 
5-invariant of the TLM. They found that 5 = 6°.3±0°.l and 5 = 6°.1±0°.2 for GX 5-1 and 
GX 17+2 respectively. But they pointed out the very large x 2 re d = 95.2/9 and xled = 377/9 
for GX 5-1 and GX 17+2 respectively. Their x 2 calculations include the statistical error 
only but do not account for the systematic resonance shift r] and for the random resonance 
error Auq. It can be concluded that their x 2 f° r the 5— mean are significantly reduced when 
this effect is taken into account. For atoll sources in general it is not clear what kind of low 
frequencies can be used for evaluation of 5. For example, vvlf presented by van Straaten et 
al. (2002) for 4U 1728-34 have much stronger dependence on v\ than those in 4U 0614+09. 
They are rather subharmonics of uvlf (than uhbo) which are interpreted by Titarchuk, 
Bradshaw & Wood (2001) as the magnetoacoustic oscillation frequencies (see also Ford & 
van der Klis 1998). 



5. Conclusions 

In this Letter, I present the consequences and predictions of the TLM as an oscillatory 
model for the QPO observations. It is fair to ask what kind of arguments one can put forth for 
the resonance effect in QPO phenomena in general. I argue that they are (i) the Lorentzian 
shape of the QPO feature (e.g. van der Klis 2000), (ii) l/oo resonance dependence of the 
rms amplitude for lower kHz QPO found by Titarchuk, Cui & Wood (2002) using Jonker's 
et al. (2000b) data; (iii) observational evidence for 2:3 resonance in the "kHz" QPO sources 
(Abramovicz & Kluzniak 2001; Abramowicz et al. 2002), (iv) harmonics and sidebands of 
low kHz QPOs (Jonker et al. 2000b). The whole issue of the observational evidence for the 
resonance effect in the QPO phenomena is beyond the scope of this short Letter. 
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Using the previously cited thorough analysis of RXTE QPO data for Sco X-l, Sco X-l, 
GX 340+0, GX 5-1, GX 17+2 and two atoll sources: 4U 1728-34, 4U 0614+09 I a detailed 
investigation of the accuracy of the QPO frequency determination in the framework of a 
weakly nonlinear oscillator is offered. The TLM as an example of a weakly nonlinear oscillator 
is used to extract the model parameters and to test the model. The TLM predicts the 
existence the invariance of the 5. It is established that: (1) The errors of the eigenfrequency 
extraction are significantly affected by the errors of the QPO halfwidth. (2) The inferred 
5— values are consistent with constants for the four Z-sources and 4U 0614+19 where kHz 
and HBO frequencies have been detected and correctly identified. I also put forth arguments 
to explain the QPO phenomena as a result of the resonance effect in anharmonic oscillatory 
systems. 

L.T. acknowledges fruitful discussions with Vladimir Krasnopolsky, Chris Shrader, Paul 
Ray and with referee. In particular, I am grateful to Charlie Bradshaw, Michiel van der 
Klis, Peter Jonker and Mariano Mendez for the data which enable me to make detailed 
comparisons of the model predictions to the observations. 
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Fig. 1. — Inferred rotational frequencies of NS magnetosphere as a function of the Kepler 
frequency for a number of Z sources: Sco X-1 (black), GX 340+0 (red), GX 5-1 (green), 
GX 17+2 (blue). Constant rotational frequencies (possibly NS spin) for Sco X-1 and GX 
17+2 (solid line) are consistent with the data. The magnetospheric rotational profiles for 
GX 340+0 (red), GX 5-1 is inconsistent with constant. The mean values of the rotational 
frequencies for these sources are given without error bars. 
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Fig. 2. — Inferred angles S between rotational configuration above the disk (possibly NS 
magnetosphere) and the normal to the plane of Keplerian oscillations as a function of the 
Kepler frequency for Z-sources: Sco X-1 (black), GX 340+0 (red), GX 5-1 (green), GX 17+2 
(blue). Constant 5— value is consistent with the data. The 5— angle (eq.[8]) is calculated 
using the lower and higher kHz peaks z/ 1; v 2 an d HBO frequencies i>hbo ( see text). 



